Fermi surface and pseudogap in electron-doped cuprates 
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Spectral functions are evaluated numerically within the t-t'-J model as relevant for electron-doped cuprates. 
The Fermi surface develops from a pocket-like into a large one with doping. The corresponding pseudogap 
in the nodal direction is well resolved at low doping and vanishes at intermediate doping. Its presence is 
strongly temperature dependent and thus connected to longer-range antiferromagnetic correlations. Features 
including double effective band and optical conductivity are consistent with experiments on electron-doped 
cuprates without invoking the closing of the Mott-Hubbard gap. 
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The insight into the interesting and challenging physics of 
high-temperature superconductors (HTSC) has been greatly 
enhanced by angle-resolved photoemission spectroscopy 
(ARPES) experiments on cuprates |1]. One of the challeng- 
ing aspects is a pronounced asymmetry of spectral properties 
between more investigated hole-doped and electron-doped 
cuprates (EDC), as revealed by recent ARPES experiments 
on representatives of the latter materials, as Nd2- 2; Ce. r Cu04 
(NCCO) HH and Sm 2 _ 2; Ce 2; Cu04 (SCCO) fl]. While the 
n-type parent compound Nd2Cu04 has very similar excita- 
tion spectra to p-type undoped cuprates, the development with 
doping is different. At low electron doping x ~ 0.04 the 
ARPES Fermi surface (FS) reveals pronounced pockets at 
k = (ir, 0) consistent with doubling of the unit cell due to 
the long-range antiferromagnetic (AFM) order 10]. With in- 
creasing x additional dispersive band-like excitation appears 
at u> < and reaches the chemical potential (u> = 0) along the 
diagonal k ~ (tt/2,tt, 2) at intermediate doping x = 0.15, 
which coincides with the disappearance of long-range AFM 
order and the onset of superconductivity. Consequently, the 
FS transforms into a large one although still with well resolved 
quasiparticles only along parts of the FS [2, 3, 4]. More recent 
analyses of NCCO y|] and SCCO at intermediate doping 
seem to resolve both effective bands and found the connection 
with the pseudogap observed in optical conductivity 
While in the latter the temperature dependence of the pseudo- 
gap reveals relation to the AFM correlations, there is so far 
no temperature-dependent ARPES data on EDC which would 
directly relate the ARPES pseudogap to AFM ordering. 

On the theoretical side, the EDC represent another very 
stringent test for microscopic models of HTSC materials. It 
has become increasingly clear that in prototype Hubbard t- 
U model and reduced t-J model additional second-neighbor 
hopping t' (as well as third-neighbor t") can account well for 
the asymmetry in several quantities |& 0]. Still, the origin 
of the FS transformation from a pocket-like to a large one as 
found by ARPES remains the central open issue. First theo- 
ries were based on the strongly renormalized interaction U 1 8] 
and the scenario of doping-decreasing U was proposed which 
results in both Hubbard bands crossing the FS at intermediate 



doping lil More recent slave-boson approach 111] and 
numerical calculations within the Hubbard model, using the 
variation-cluster |12] and dynamical-cluster approaches 11311 
give indication that most features of FS reconstruction can 
be reproduced with doping-independent U. More detailed 
insight into the spectral properties has been obtained from 
the numerical exact-diagonalization (ED) study of t-t'-t"-J 
model at T = |0]. However, results show persistence of 
the pseudogap in the nodal direction even in the overdoped 
regime. Also, the temperature behavior of the pseudogap in 
spectral function (SF) has not been clarified yet. 

In this Letter we present results of a numerical study of the 
t-t'-J model at finite T > which reveal that most qualita- 
tive features of ARPES results are a nontrivial consequence 
of strong correlations and do not require a doping-dependent 
model parameters, in particular not the closing of the Mott- 
Hubbard gap. We show that the pseudogap vanishes with 
doping in the nodal direction as well as above the pseudo- 
gap temperature T* being closely tied to the disappearance of 
longer-range AFM correlations. In the same way the FS trans- 
forms from a pocket-like into a continuous large one. Further, 
our results allow for the interpretation of two effective sub- 
bands both emerging from an AFM splitting of a correlated 
(upper Hubbard) band. Numerically evaluated a(oj) can be 
also explained within the same pseudogap scenario. 

In the following we study the t-t'-J model 
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where cj are projected fermionic operators not allowing the 
double occupancy of sites. On a square lattice we include 
besides the nearest-neighbor hopping f y = t also the second- 
neighbor hopping tij = t', where we take t' = 0.3t, as rele- 
vant for EDC • We also fix J = 0.3 t. 

For the calculation of the SF, A(k,uj), we use the finite- 
temperature Lanczos method (FTLM) 11411 . It has been shown 
that on small systems reachable by ED the calculation at 
T > gives more macroscopic-like results and gives a better 
controlled approach to T — > provided that T > Tf s where 
Tf s is the finite-size temperature determined by the model and 
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system size |[l4ll . For cases considered here with tilted square 
lattices of N = 18, 20 sites and iV e = — 4 doped electrons 
we have TV S ~ 0.1 t. On small systems with fixed boundary 
conditions one can consider only a discrete set of wavevec- 
tors k;, I = l,N. To scan the whole Brillouin zone of k we 
employ twisted boundary conditions by introducing hopping 
elements — > ty = tij cxp(i8 ■ fy) in Eq. (Q]), in this way 
reaching arbitrary momenta k = k; + 9 [7]. Using FTLM we 
calculate the Green's function G(k, u>) and therefrom extract 
the self energy S(k, u>) = uj — (k — a/G(\c,u>). Borrow- 
ing the idea from cluster dynamical mean-field approaches we 
smoothen £(k, u>) over small interval 5k ~ 0.3 and then re- 
calculate SF. For details we refer to Ref.l 161. 




Figure 1: (Color online) Contour plot of A(\i,u) — 0) representing 
the Fermi surface for different electron doping: (a) c e = 1 /20, (b) 
c e = 2/20, (c) c e = 3/20, and (d) c e = 4/18. Results are shown 
for T/i = 0.1. 



Let us start with the presentation of results for the FS evo- 
lution. In Fig. 1 we display the scans of the SF at the chemical 
potential ^4(k, oj = 0), as usual also in ARPES presentations 
of FS HI 0, IH 0], in the quarter of the first Brillouin zone for 
different electron dopings c e = N e /N. Data in Fig. 1 are sup- 
plemented with Fig. 2 where full A(k, u>) are presented along 
the AFM zone-boundary line. 

At the lowest doping c e = 1/20 in Fig. la we get quasipar- 
ticle weight forming a circle-like structure around the (ir, 0) 
point, well known as the electron pocket consistent with sev- 
eral experimental and theoretical studies JH 0, S U]. The 
FS scan looks similar also for c e = 2/20 in Fig. lb except 
that the pocket gets more asymmetric weight (opening from 



In this doping regime the FS can be defined in a usual way, i.e. 
as a k-line of points where the quasiparticle peak in SF crosses 
the chemical potential, represented as white dots around (tt, 0) 
in Figs. la,b. On the other hand, white crosses in Figs, la-c 
indicate a k-line where the SF shows symmetric pseudogap. 
Note that the line of zeroes in the real part of the Green's func- 
tion determining the Luttinger sum rule is then represented by 
the combination of both lines. It is indicative in Fig. lb that 
already for c e = 2/20 both lines nearly touch as a precursor 
of a large FS formation. 

As the doping is further increased, in Fig. lc, the FS is 
formed almost along the whole Luttinger line. The pockets 
are completely opened and the transfer of the SF weight into 
the nodal region is observed, although the pseudogap feature 
remains visible as seen in Fig. 2c. Finally, at c e — 4/18 
in Fig. Id the pseudogap disappears also in the nodal direc- 
tion resulting in a large continuous FS. In addition, the SF in 
Fig. 2d are well consistent with the usual quasiparticle peak 
crossing the chemical potential. 




the k y — side) around (jr, 0). Within the rest of the Bril- 
louin zone, in particular along the AFM zone-boundary, the 
SF shows a pseudogap at uj ~ 0, as evident also in Figs. 2a,b. 



Figure 2: Spectral functions A(k, uj) along the AFM zone-boundary 
line. Dopings are the same as in Fig. 1 and T/t = 0.1. 

In Fig. 3 we focus in more detail on the doping and tem- 
perature dependence of the pseudogap. We present A(k, io) 
in the nodal point of the Luttinger surface located as shown in 
the inset of Fig. 3a. It is evident that a pronounced dip in SF 
at low doping c e < 0.15 fills up with increasing doping and 
transforms into a usual peak at c e = 4/18. The pseudogap 
is clearly manifested also in the behavior of the correspond- 
ing self-energies S"(k, u>) in Fig. 3b. Besides the overall lin- 
ear dependence | E" (k, oj) | oc li for uj > 0, characteristic for 
anomalous (marginal) Fermi liquid and well established also 
in hole-doped cuprates HI [l6Tl . additional narrow Lorentzian 
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at ui ~ reveals the pseudogap and its strength fiei a gain, 
this feature vanishes at c e = 4/18, but as well disappears 
with increasing T. The intensity (the area) A 2 of the pseudo- 
gap contribution in S"(k, ui) is the measure of the pseudogap 
width, i.e. oc A II 1611 . Therefore we present in Fig. 3c T- 
dependence of A 2 for different c e . It is evident that temper- 
ature T* above which the pseudogap effectively disappears 
decreases with c e and exists only in the underdoped regime 
with c e < 4/18. This clearly links pseudogap to the presence 
of longer-range AFM correlations in the system. In this sense, 
the pseudogap discussed here is related to large (high-energy) 
pseudogap observed in hole-doped cuprates yj]. 
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Figure 3: (Color online) a) A(k,u>) in the nodal point of the Lut- 
tinger surface presented in the inset, b) corresponding self-energies 
E"(k, ui), and c) pseudogap intensity A 2 vs. T for different c e . 

In Fig. 4 we present a weight map of A(k, ui) along sym- 
metry lines in the first Brillouin zone. Again, one can notice 
strong electron-pocket contribution close to (w, 0) and evident 
pseudogap in the nodal direction at lowest c e = 1/20 which 
becomes less pronounced at higher doping and disappears at 
c e = 4/18. In addition, it is interesting to follow effective 
bands which are clearly seen in Fig. 4. At lowest c e = 1/20 
one can well resolve the two-band structure, at least along cer- 
tain k directions, e.g., along (0, 0) - (ir/2, tt/2), (0, 0) - (tt, 0) 
as well as (n, 0) — (n/2, tt/2). This feature dissolves with 
doping and finally at c e = 4/18 transforms into a single band. 
Such a two-band structure is very close to the one observed 
and analyzed in ARPES on EDC The weaker band can 

be interpreted as the shadow-one due to AFM order and its 
intensity can be linked to AFM order parameter s. In fact, re- 
sults in Fig. 4 can be well represented by effective bands of 
the form, 



where 7k = (cosfc x + cosfc y )/2, 7^ = cos k x cos k y are 
corresponding to first and second-nearest neighbor effective 
hopping t, t', respectively. The dispersion of the form © can 
be considered as an effective tight-binding band split due to 
long-range AFM s 7^ and can be analytically derived from 
the t-t'-J (or Hubbard) model using the memory-function ap- 
proach 1 15] whereby t, t' are renormalized by strong correla- 
tions. Very similar bands have been recently used to describe 
observed ARPES spectra in SCCO fl. 

In Fig. 4 we also plot fits of e±(k), with corresponding 
intensity of subbands shown in the scale from black (no in- 
tensity) to white. We note that the main doping dependence 
emerges through s — 0.35 — which is the largest for low 
doping and zero at c e = 4/18 while other parameters are only 
weakly doping dependent, e.g. t/t = 0.3 — 0.6,— t'/t = 
0.1 -0.2 and w = t 2 . 
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Figure 4: (Color online) Weight map A(k, ui) along symmetry lines 
in the first Brillouin zone for different dopings and T/t = 0.1. The 
fit to split subbands is shown with the corresponding intensities. 

Here we comment also on the asymmetry of high-energy 
dispersion (u> < 0) between hole and electron-doped sys- 
tems. From Fig. 4 we can clearly see that the dispersion from 
(7r/2, 7t/2) to (0, 0) is a rather standard one without signifi- 
cant kinks, but for u> > 0, which corresponds to u> < case of 
the hole-doped system in t — t' — J models, a kink behavior 
is seen along (7r/2, 7r/2) to (7r, 7r) as well as from (tt,0) to 
(7r, 7r). This difference is consistent with related asymmetry 
observed in ARPES recently lfl8ll . The details will be given 
elsewhere 1 19]. 

Finally, let us discuss the relation with the optical conduc- 
tivity a{ui) in EDC ||4j,|5J which also shows pronounced pseu- 
dogap at low doping. First, we evaluate cr(ui) directly by 
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calculating the dynamical current-current correlations in the 
ground state at T = and via FTLM at T > QQ3]. Re- 
sults for T = and T/t = 0.15 t at different c e are presented 
in Fig. 5. On the other hand, to establish the relation with the 
pseudogap observed in SF we also consider the simplest (de- 
coupling) approximation in terms of SF neglecting the vertex 
corrections, 

o-dM = ^E(^) 2 / du'g(cj,w')A(k,u')A(k,u)'-u), 

k 

(3) 

where g(u),u>') — f(uo' — u>) — /(w') and Vk are bare band 
velocities as determined by t, t'. Numerically evaluated SF as 
presented above are inserted into Eq. <£3j and results are plot- 
ted in Fig. 5 as black lines. A nontrivial message following 
from Fig. 5 is that <7d(w) yields very good overall qualitative 
and even quantitative agreement with the exactly evaluated 
<t(oj) (at T > 0), whereby features are plausibly broader in 

CTd(w). 
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Figure 5: (Color online) Optical conductivity <j(uj) as obtained from 
the small-system calculation at T = (red line), shown for c e = 
1/20, 2/20, 3/20, and T/t = 0.15 (blue line). Black lines represent 
Od{u) calculated from the convolution of spectral functions at T/t = 
0.1. 

At low doping, Ch = 1/20, 2/20 we can observe at T/t = 
0.15 t a hump in a(to) at uj/t ~ 0.7, even more pronounced 
at T = 0, which can be interpreted as a manifestation of 
the pseudogap. The same peak-dip feature is also resolved 
(although with somewhat shifted peak) in convoluted ad{u>) 
which is the indication that the pseudogap in conductivity has 
the same origin as the one in the SF. Furthermore, vu is the 
largest along the nodal part of the FS so it is plausible that 
<j(lo) also tests the nodal pseudogap. Further doping sup- 
presses the peak-dip feature in Fig. 5 at c e = 4/18 which 
can be related to the simultaneous pseudogap disappearance 
in A(k,uj). 

In conclusion, presented results for spectral functions and 
optical conductivity show overall good qualitative and even 
quantitative agreement with experimental results on EDC with 



different doping, both for ARPES on NCCO and SCCO 
|E El 01, and cr(w) JH Etl- Since we have used the simplest 
prototype model, Eq.©, for strongly correlated electrons with 
fixed parameters, observed nontrivial phenomena as the Fermi 
surface reconstruction, pseudogap and split bands appear as a 
direct consequence of strong correlations. So there is no need 
to invoke a doping-dependent Mott-Hubbard gap. 

In comparison with more investigated hole-doped cuprates, 
EDC seem to have an advantage that anomalous SF as well as 
<j(lo) can be explained in a more straightforward way by in- 
voking the long-range AFM order s, persisting in EDC mate- 
rials up to c e = 0.15. In particular, the double effective bands 
as seen in Fig. 4 at lower doping as well as in recent ARPES 
experiments lil 01 can be interpreted with a band splitting 
proportional to s, Eq.©, which disappears in the overdoped 
regime. 

The pseudogap which is pronounced in SF along the zone 
diagonal and in <t(uj), shows strong doping and temperature 
dependence, whereby we can establish the pseudogap temper- 
ature T*(c e ) consistent with the one deduced from the optical 
conductivity a (10) (5J] but so far not analyzed with ARPES. 
Both doping and T-dependence clearly establish longer-range 
AFM correlations as the origin of the pseudogap in EDC. 
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